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Abstract 

Consider real symmetric, complex Hermitian Toeplitz and real sym- 
metric Hankel band matrix models, where the bandwidth b]\[ ~* oo 
but bisi/N ^ b, b G [0, 1] as iV ^ oo. We prove that the distributions 
of eigenvalues converge weakly to universal, symmetric distributions 
7^(6) and ^„{b). In the case 6 > or 6 = but with the addition of 
&Af > C N2+'^° for some positive constants eo and C, we prove almost 
sure convergence. The even moments of these distributions are the 
sum of some integrals related to certain pair partitions. In particular, 
when the bandwidth grows slowly, i.e. 6 = 0, 7r(0) is the standard 
Gaussian distribution and 7^(0) is the distribution \x\ exp(— a;^). In 
addition, from the fourth moments we know that the 7^,(6) 's are dif- 
ferent for different 6's, the 7^(6)'s different for different b E [0, i] and 
the 7^(6)'s different for different b G [5, 1]. 



1 Introduction 



In Random Matrix Theory, the most important information is contained 
in the eigenvalues of matrices and the most prominent analytical object is 
the distribution of eigenvalues. That is, for a real symmetric or complex 
Hermitian N x N matrix A with eigenvalues Ai, • • • ,Xn, the distribution of 

°Key words: Random matrix theory; distribution of eigenvalues; Toeplitz band ma- 
trices; Hankel band matrices. 
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its eigenvalues is the normalized probability measure 

1 ^ 

^'^■■=NE^^r (1-1) 

If ^ is a random matrix, then fi^ is a random measure. 

The limit distribution of (jl.ip is of much interest in Random Matrix 
Theory. In [18j, Wigner found his famous semicircular law and then in 
pointed out that the semicircular law was valid for a wide class of real sym- 
metric random matrices. Since then much more has been done on various 
random matrix models, a standard reference is Metha's book [IB]. Recently 
in his review paper [1], Bai proposes the study of random matrix models 
with certain additional linear structure. In particular, the properties of the 
distributions of eigenvalues for random Hankel and Toeplitz matrices with 
independent entries are listed among the unsolved random matrix problems 
posed in [1], Section 6. Bryc, Dembo and Jiang ^ proved the existence 
of limit distributions 7^ and 7^, for real symmetric Hankel and Toeplitz 
matrices. The moments of 7^ and 7^, are the sum of volumes of solids, 
from which we can see that 7^ and 7^ are symmetric and of unbounded 
support. At the same time Hammond and Miller [12] also independently 
proved the existence of limit distribution for symmetric Toeplitz matrices. 
In the present paper we shall prove the existence of limit distribution of 
eigenvalues for real symmetric and complex Hermitian Toeplitz band ma- 
trices, and also for real symmetric Hankel band matrices. Random band 
matrices have arisen in connection with the theory of quantum chaos [H [11] 
, and the limit distribution of eigenvalues of band matrices has been studied 
in [31 [IT]. 

Note that Toeplitz matrices emerge in many aspects of mathematics 
and physics and also in plenty of applications, e.g. p], especially [21 [9] for 
connections with random matrices. Hankel matrices arise naturally in prob- 
lems involving power moments, and are closely related to Toeplitz matri- 
ces. Explicitly, for a Toeplitz matrix of the form T/v = and 
a Hankel matrix of the form Hn = (aj+j„2)i^=i5 if Pn = {^i-i,N-j)i^j=i 
the "backward identity" permutation, then PnTm is a Hankel matrix and 
Pj\iHn is a Toeplitz matrix. In this paper we always write a Hankel matrix 
Hj^ = PnTn where we assume the matrix entries a_Ar+i, • • • , ao, • • • , a^-i 
are real- valued, thus Hf^f is a real symmetric matrix. In addition, if we intro- 
duce the Toeplitz or Jordan matrices B = and F = {dij+i)fj^i, 
called the "backward shift" and "forward shift" because of their effect on the 
elements of the standard basis {ei, • • • , cat} of M^, then an N x N matrix 
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A can be written in the form 



N-l AT-l 

A=Y^ a-jB^ + J2 flj^^ (1-2) 

3=0 j=l 

if and only if A is a Toeplitz matrix where a_Ar+i,-- - ,ao,--- ,aAf-i are 
complex numbers |13j . It is worth emphasizing that this representation of a 
Toeplitz matrix is of vital importance and the starting point of our method. 
The "shift" matrices B and F exactly present the information of the traces. 
Consider a Toeplitz band matrix as follows. Given a band width < N, 

let 

r] - = 1"^' - (1 3) 

1 0, otherwise. 

Then a Toeplitz band matrix is 

Ttv = aj_j)fj=i. (1.4) 
Moreover, the Toeplitz band matrix Tjv can be also written in the form 

Tn = Y^ a^jB^ + ajFK (1.5) 
i=o j=i 

Obviously, a Toeplitz matrix can be considered as a band matrix with the 
bandwidth = N — 1. Note that when referring to a Hankel band matrix 
Hn, we always mean Hf^ = PnTn where Tf^ is a Toeplitz band matrix. 
In this paper the three models under consideration are : 
(i) Hermitian Toeplitz band matrices. The model consists of A'^-dimensional 
random Hermitian matrices T^q = {r]ij ai-j)^j_^ in Eq. (jl.4p . We assume 
that Reoj = Rea_j and Imaj = — Imo_j for j = 1,2, and {ao} U 
{ReOj,Imaj}jgN is a sequence of independent real random variables such 
that 

E[aj] = 0, E[|ajf ] = 1 for j eZ (1.6) 



and further 



supE[|ajn = Cfc < oo for G N. (1.7) 



Notice that 



supE[|aj|''] < oo sup{E[|Reaj|''],E[|Imaj|'']} < oo. 
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and the assumption (|1.7p shows that when fixing k, these moments of the 
independent random variables whose orders are not larger than k can be 
controlled by some constant only depending on k. 

The case of real symmetric Toeplitz band matrices is very similar to the 
Hermitian case, except that we now consider A'^-dimensional real symmetric 
matrices Tn = (% aj-i)y=i • 

(ii) Symmetric Toeplitz band matrices. We assume that aj = a-j for 
j = 0,1, ■ ■ ■ , and {aj}j^Q is a sequence of independent real random variables 
such that 

E[aj]=0, E[|ajf ] = 1 for j = 0,l,---, (1.8) 

and further 

supE[|aj|''] = Cfe < oo for G N. (1.9) 

i>o 

(iii) Symmetric Hankel band matrices. Let = {rjij ai-j)^^^-^ be a 
Toeplitz band matrix and = PnTn be the corresponding Hankel band 
matrix. We assume that {aj}j^z is a sequence of independent real random 
variables such that 

E[aj]=0, E[|ajp] = l for jeZ (1.10) 

and further 

supE[|aj|''] = Cfc < cx) for A; G N. (1.11) 

We assume that bj\f grows either as 1) b]\f/N — > 6 as iV ^ oo, 6 G (0, 1] 
(proportional growth), or as 2) ^ cxd as ^ oo, b]\f = o{N) (slow 
growth). Under these assumptions we will establish the limit distributions 
7j,(6) for Toeplitz band matrices and ^^{h) for Hankel band matrices. In the 
case 6 > 0, we prove that the distribution of eigenvalues for Toeplitz (Han- 
kel) band matrices converges almost surely to 7-j,(6) {'^^{h)), the moments 
of which are the sum of some integrals depending on h. For 6 = 0, we obtain 
that 7j, (0) is the standard Gaussian distribution and 7^ (0) is the distribu- 
tion |x| exp(— x^). In the case 6 = 0, with the addition oib^s >C N^'^^'' for 
some positive constants eo and C, we also prove almost sure convergence. 

The plan of the remaining part of our paper is the following. Sections 
2 and 3 are devoted to the statements and proofs of the main theorems for 
band matrices with proportionally growing and slowly growing bandwidths 
b^, respectively. In Section 4, the fourth moments of the limit distributions 
7j,(6) and 7^(6) are calculated, which shows their difference for different 6's. 
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2 Proportional Growth of 67V 



In order to calculate the moments of the limit distribution, we review some 
basic combinatorical concepts. 

Definition 2.1. Let the set [n] = {1, 2, • • • , n}. 

(1) We call vr = {Vi, • • • , Vr} a partition of [n] if the blocks Vj {I < j < r) 
are pairwise disjoint, non-empty subsets of [n] such that [n] = Vi U • • • U V^. 
The number of blocks of vr is denoted by |7r| , and the number of the elements 
of Vj is denoted by | V^- 1 . 

(2) Without loss of generalization, we assume that Vi, • • • ,Vr have been 
arranged such that si < S2 < ■ ■ ■ < Sr, where Sj is the smallest number of 
Vj. Therefore we can define the projection 7r(i) = j if z belongs to the block 
Vj] furthermore for two elements p, q of [n] we write p q if T^ip) = 7r(g). 

(3) The set of all partitions of [n] is denoted by V{n), and the subset 
consisting of all pair partitions, i.e. all \Vj\ = 2, 1 < j < r, is denoted 
by V2{n). The subset of 'P2{'n) consisting of such pair partitions that each 
contains exactly one even number and one odd number is denoted by V\{n). 
Note that V2{n) is an empty set if n is odd. 

We can formulate our results for Hankel and Toeplitz band matrices with 
the proportional growth of as follows. By convention, Ib represents the 
characteristic function of a set B. 

Theorem 2.2. Let he either a Hermitian ( U.6\) - (1.7^ ) or real symmet- 
ric ( lil.8\ )- [1.9\) ) Toeplitz random band matrix, where bj^/N b as N ^ 00, 
b £ (0,1]. Take the normalization = T/v/^/ (2 — b)bN , then fix^ con- 
verges almost surely to a symmetric probability distribution Jj,{b) which is 
determined by its even moments 

^ ^ 7re-P2{2fc) ^"'-^J^' ^'-^J j=l i=l ;=0 

(2.1) 

1 if i is the smallest number of tt ^(7r(f)); otherwise, €7^(1) = 

Theorem 2.3. Let Hj\f be a real symmetric ( lll.l(J\) - lll.ll\) ) Hankel random 
band matrix, where h^ /N b as N ^ 00, b £ (0, 1]. Take the normalization 
Y]\f = Hj\f/y^ (2 — b)bN, then hyj^ converges almost surely to a symmetric 



m2k{lT{b)) = 

where €7^(2) = 
-1. 
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probability distribution (b) which is determined by its even moments 

2k j k 



(2-^)\,^,,)^[0,i]x[-i,i]%i-\ to 

(2.2) 

Let us first give two basic lemmas about traces of Toeplitz and Hankel 
band matrices. Altliough tlieir proofs are simple, they are very useful in 
treating random matrix models closely related to Toeplitz matrices. 

Lemma 2.4. For a Toeplitz band matrix T = {rjij ai^j)fj^i with the hand- 
width bM where a-N+i,-" ,ciN-i o,re complex numbers, we have the trace 
formula 

N bN k k I 

i=l ji,- ,jk=-bN 1=1- 1=1 ?=1 ^'^1^" 

(2.3) 

Proof. For the standard basis {ei, • • • , e^} of the Euclidean space M^, we 
have 

6]v fejv fejv 

Tei = '^a-jB^ Ci + '^ajF^ ei= ^ ajl^i^^^ii + j) Ci+j. 
j=o i=i j=-bN 

Repeating T's effect on the basis, we have 

bjv k k I 

ji,- ,3k=-bN 1=1 1=1 q=l ''^^1^" 

k 

By tr(T'=) = e\T^ e,, we complete the proof. □ 



i=l 



Lemma 2.5. Given a Toeplitz band matrix T = {r]ij ''H-j)ij=i with the 
bandwidth b^ where a_iv+i, • • • ,aiv-i are real numbers, let H = PT he the 
corresponding Hankel band matrix where P = {5i-\^N-i)f^j=i- We have the 
trace formula 
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tr(T^) = 

N bpf k k I 

E E n flj; n %,Af](* - E ) A: even, 

i=i ii,- jfc=-bjv«=i «=i 9=1 0, E(-i)«i9 

Af 6]v fc A: / 

E E na^, n%,^](^- E(-i)%) , kodd. 

(2.4) 

Proof. Follow a similar procedure as in the proof of Lemma 12.41 □ 

We are now ready to prove the main results of this section. 

Proof of Theorem \2.SX We prove the theorem by the following steps: 
Step 1. Calculation of the Moments 
Let 

mk,N = nj x^t^xj^idx)] = lE[tr(X^)]. 
Using Lemma 2.4, we have 

N bN k I k 

mk,N= Yi n%.^i(^+z]-?9) *^ . ^[^"J'^ 
For j=(ji, • • • ,jk), we construct a set of numbers 5j = {|ji|, . . . , \jk\} with 
multiplicities. Note that the random variables whose subscripts have differ- 
ent absolute values are independent. If 5j has one number with multiplicity 
1, by independence of the random variables, we have ]E[]^j!^^ aj^] = 0. Thus 
the only contribution to the k-th moment comes when each of Sj has mul- 
tiplicity at least 2, which implies that Sj has at most [|] distinct numbers. 
Once we have specified the distinct numbers of 5j, the subscripts ji, • • • ,jk 

tions (fLT]) and (fL9]) . we find 



are determined in at most 2'^([|])'^ ways. By independence and the assump- 



mk,N = 0((6jv)"^+'^') = 0(iV~t+[|l 



Therefore, for odd k 



lim mfc AT = 0. 

N-*oo 



It suffices to deal with m2fe, AT. If each of Sj = {| ji|, • • • , |j2fc|} (here j=(ji, • • • ,j2fc)) 
has multiplicity at least 2 but one of which at least 3, then Sj has at most 
k — 1 distinct numbers. Thus, the contribution of such terms to m2k n is 
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0{N^^). So it suffices to consider all pair partitions of [2A;] = {1, 2, • • • , 2k}. 
That is, for vr G V2{2k), if p q, then it is always the case that jp = jq or 
jp = —jq. Under the condition 

2k 
9=1 

according to (j2.3p . considering the main contribution to the trace, we should 
take jp = —jq ; otherwise, there exists po,qo S [2k] such that 



Jpo ~ 3qo 2 



^ 2k 



We can choose other A: — 1 distinct numbers, which determine = jq^. 
This shows that there is a loss of at least one degree of freedom and the 
contribution of such terms is 0{N^^). Thus aj^ and aj^ are conjugate. So 
we can write 

j^-k-l N bN 2k I 

^ ' ■K&V2{2k) i=l h,- ,jk=-bN 1=1 9=1 

(2.5) 



Follow a similar argument from [71[3]: for fixed vr E 7^2 (2fc), 

Y \{h,N]{i + Y^-^'i^M<i))^ 

i='^ jl,--- ,jk=-bN 1=1 9=1 

i.e. 

N bN 2k . I 

«=1 Jir-- Jfc=-f'iv ^=1 9=1 

can be considered as a Riemann sum of the definite integral 

/ n^[o.i](^o +xi^'^(^)^'r(9))n^^'- 

i[o,i]x[-M'=fi t^i to 

Given vr G V2{2k), as N ^ oo, each term in ()2.5p can be treated as an 
integral. Thus we have 7712^(7^(6)) = lim m2k,N is the representation of 

A^— >oo ' 



the form (pTj) . 

Step 2. Carleman's Condition 
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Obviously, from ()2.ip we have 

^ k 

[2-bf ^±^,,,7[o,iixhi,il* (2-6)'= 



TT£_V2{2k) 



By Carleman's theorem (see (lOj), the hmit distribution 7^(6) is uniquely 
determined by the moments. 

Step 3. Almost Sure Convergence 

It suffices to show that 

oo 

^E[(tr(X^) - E[tr(X^)])^] < oc (2.6) 

N=l 

for each fixed k. Indeed, by Lemma 12.41 we have 

tr(X^) = ^ J^Abo;?] 

N2 

PO,P 



^[po; p] = n n hhN] (po + X] Pi")^ 

1=1 1=1 q=l 



where 

k k I 

" fe 

E 

9 = 1 

p = (pi, . . . ,pfc), and the summation p runs over all possibilities that 
p G {-hN, . . . , hN}^ and po G {1, • • • , A^}- Thus, 

^E[(tr(X^)-E[tr(X^)])4] 

= ^ ~J,^2k — E n i^po-M-mWM])] (2.7) 

»o.io.so.*o pe{i,j,s,t} 

where the summation X]'o>io.so>'o runs over all possibilities that i,j,s,t G 

{-bN, ■ ■ .,hNY and icJcSc^o G {1, • • • ,^}- 

For i = (ii, . . . , ifc), as in Step 1 we still construct a set of numbers 
S\ = {Nil) • • • ) Kfc|} with multiplicities. Obviously, if one of 5i, . . . , ^t, for 
example, Si does not have any number coincident with numbers of the other 
three sets, then the term in (j2.7p equals by independence. Also, if the 
union S = S\IJ ■ ■ - VJ St has one number with multiplicity 1, the term in ()2.7p 
equals 0. 

Now, let us estimate the non-zero term in ()2.7p . Assume that S has 
p distinct numbers with multiplicities i/i, . . . ,fp, subject to the constraint 
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ui + . . . + Up = 4k. If p < 2A; - 2, by ([L6])-(lL9]), the contribution of the 
terms corresponding to S is bounded by 

(2 _ ^)-2fcA-2fc 



for some constants and Ckfi- 

For the case of p = 2A; — 1, we have two numbers shared by three of four 
sets Si, . . . ,St each or one number shared by the four sets. It is not hard 
to see that there always exists one of Si, . . . , St, for example, Si, which has 
one number with multiplicity 1, denoted by for some go G {1) ■ • ■ > 
Consequently, under the constraint 

k 
9=1 

k 

according to (j2.3p . we have ig^ = iqo — Hence we can choose the other 

9=1 

2k — 2 distinct numbers of S, which determine the subscript ig^. This shows 
that there is a loss of at least one degree of freedom and the contribution of 
the terms corresponding to S is bounded by CkfiN~'^. 

The case oi p = 2k implies each number in S occurs exactly two times. 
Thus there exist two of Si,...,St, for example. Si and 5j, which share 
one number with one of S^ and St respectively. Consequently, under the 
constraints 

k k 

^ = and ^ = 

q=l q=l 

according to ()2.3p . there is a loss of at least two degrees of freedom and the 
contribution of the terms corresponding to S is bounded by Ck^b^~'^- 
Therefore, 

^E[(tr(X^) - E[tr(X^)])^] < Ck,tN-' 

and almost sure convergence is proved. Consequently, the proof of Theorem 
12.21 is complete. □ 

Proof of Theorem \2.3[ . Using Lemma 12.51 and our assumptions, the proof is 
very similar to that of of Theorem 12.21 Here we don't repeat the process. 
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Note that in the representation of even moments in Eq. ()2.2p the pair parti- 
tion vr G 7^2(2^) comes from the constraint 

k 

E(-1)'J. = (2.8) 

q=l 

in the form of the trace for Hankel matrices in Lemma 12.51 □ 

3 Slow Growth of 

We give our results for Hankel and Toeplitz band matrices with slow growth 
of the bandwidth as follows. 

Theorem 3.1. LetT]\j he either a Hermitian ( il.6\) - [1.7^ ) or real symmetric 
( il.8\) - D.y\) ) Toeplitz random band matrix, where bjv — > 00 but b]\f/N 
as N ^ 00. Take the normalization = T^/^/Tb^, then ^Xpf converges 
weakly to the standard Gaussian distribution. In addition, if there exist some 
positive constants eo and C such that b^ > CA^a+'^o^ then converges 
almost surely to the standard Gaussian distribution. 

Theorem 3.2. Let H^^ be a real symmetric ( il.lO\) - il.ll\) ) Hankel ran- 
dom band matrix, where b^ ^ 00 but bj^/N —> as —> 00. Take the 
normalization = H]\f 6jv, then /iy^^ converges weakly to the distribu- 
tion f{x) = |x|exp(— x^). In addition, if there exist positive constants eo 
and C such that b^ > C Ni^'^° , then converges almost surely to the 
distribution f{x) = \x\ exp(— x^). 

Remark 3.3. For real symmetric palindromic Toeplitz Matrices, i.e. real 
symmetric Toeplitz matrices under extra conditions: aj_i = a7v-j,0 < j < 
N, Massey, Miller and Sinsheimer [T5] have obtained the Gaussian normal 
distribution for eigenvalues. And for random reverse circulant matrices, i.e., 
Hankel matrices under extra conditions: a_j = aAr_j, 1 < j < A^, Bose and 
Mitra [5j obtained the same distribution of eigenvalues f{x) = \x\ exp(— x^) 
as in Theorem 13.21 

Since proofs of both theorems are very similar, we prove only the Toeplitz 
case. 

Proof of Theorem \3.1[ The proof is quite similar to that of Theorem 1 2 . 2 1 and 
some details will be omitted. We will lay a strong emphasis on the derivation 
of the Gaussian distribution. Here the slow growth of bN leads to an easy 
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calculation of the complicated integrals. We now complete the proof of the 
Gaussian law by showing 

(1) mk^N = ;^lE[tr(X^)] converges to the k-th moment of the standard 
Gaussian distribution. 

(2) Assume that > C A^a+'^o for some positive constants eo and C. 
For each fixed /c, 

^E[(tr(X^) - E[tr(A^)])^] < oo. (3.1) 

N=l 

By Lemma [23] and our assumptions ()1.6p - ()1.9p . it follows that 

mfc,;v = 0((6^)-t+[tl). 
Since b]\f ^ oo as N ^ oo, for odd k 

lim mfc AT = 0. 

N—i-oo 

It suffices to deal with m2k,N- However, again by Lemma 12.41 and our as- 
sumptions, the contribution with the exception of all pair partitions to m2k,N 
is 0{b~^^). So it suffices to consider all pair partitions of [2A;] = {1, 2, • • • , 2k}. 
That is, for tt G V2{2k), if p q, then it is always the case that jp = jq or 
jp = —jq. Under the condition 

2k 
q=l 

according to (12. 3p . considering the main contribution to the trace, we should 
take jp = —jq (otherwise, there is a loss of at least one degree of freedom). 
Thus and aj^ are conjugate. So we can write 

N bN 2k I 

^ 7re-P2(2fc) i=l ji,--- ,jk=-bN «=1 9=1 

(3.2) 

Note that the constraints 

1 <i + J^e.(g)i.(g) < A, l</<2/c, 

q=l 

i.e. ^ 

-<^ + ye,(g):^<l, 1< l< 2k. 

9=1 
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Since 

I 



I / .fcTrWj 
9=1 

as N ^ oo, we have 



TTvy; ^ I _ ^ 

q=l 



I 

h,N\{i + ^(^n{q)jn{q)) = 1 + 0(1). 

9=1 

Hence 

m2k,N = 0(1) + {2k - 1)!! (^^g±l)^(l + o(l))2fc (3.3) 

converges to {2k — 1)!! as ^ 00. Assertion (1) is then proved. 

We now prove (2). The same argument as in the proof of almost sure 
convergence in Theorem 12.21 shows that 



ijE[(tr(X^)-E[tr(X^)])4] = 0(^; 



Since bj^ > C A^a+'^o, we have 

^E[(tr(X^) - E[tr(X^)])4] < C,^^, 

where is a constant depending on k only. This completes the proof of 
assertion (2). 

The proof of Theorem 13. II is then complete. □ 



4 First Four Moments of ^j,{b) and 7^(6) 

We will compute the second and fourth moments of 7j,(6) and ^^{h). From 
the fourth moment we can read that the 7j,(6)'s are different for different 
6's, the 7^(6)'s different for different h e [0, ^] and the 7^(6)'s different for 
different h € [^,1]. However, for the 2/c-th moments {k > 3), the integrals 
on the right-hand sides of (j2.ip and (j2.2p are in general different for different 
pair partitions, which makes difficult the explicit calculations of the higher 
moments. 

Observe that for vr E 'P\{2k), the corresponding integrals on the right- 
hand sides of ()2.ip and (j2.2p are in fact the same. Therefore, it is sufficient 
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to calculate the integral in the Toeplitz case for every pair partition. For 
the pair partition vr G V2{'^k), we introduce the symbol 



.2k j k 

PTvib) = / ]T/ro,i](xo + 6Ve^(i)x^(j))TTdxz. 



For k=l it is easy to obtain the second moments 

m2{jj.{b)) = 1, 7712(7^(6)) = !. 

Thus the fourth moment is the first "free" moment in seeing the shape of the 
distribution. When A; = 2, a direct calculation (this can be checked using 
Mathematica) for all pair partitions 

vri = {{l,2},{3,4}},7r2 = {{1, 4}, {2, 3}}, ^3 = {{1, 3}, {2, 4}} 

of 7^2 (4) shows 

JiW (4.1) 



-l+6fe--2fe3 1,^(1 11 



and 



P^sW - \ 2(-l+6fe-6fe2+263) ,1 l4.2j 

I 3P ' Ofcl2'J^J- 

Thus, for b G [0, i], 



'^4(7t(o)) ~ 



3(2 - 6)2 

which strictly decreases on [0, ^]. When b G (^, 1], 

4(-l + 66-362) 



m4(7T(^)) 



362(2 - 6)2 



which also strictly decreases on (^,1]. 

Therefore, one knows that 77^4(7^(6)) strictly decreases on [0, 1]. Further, 
the distributions 7j,(6)'s are different. In particular, 7j,(6) (0 < 6 < 1) is 
indeed different from the normal distribution. Note that for 6=1 the fact 
that the limit distribution is not Gaussian has been observed in [41 [71 [12] . 

We turn to the Hankel-type distribution 'Jj^ib). From ()4.ip . one obtains 
that for 6 G [0, ^] 

4(6 - 56) 
"^4(7.(6)) = ^^^3^ 



14 



which strictly decreases on [0, ^] while for 6 G 1], 

, ,,,, 2(-l + 66-262 

"l4(7ff(0)) = 



362(2-6)2 

which strictly increases on (^,1]. Thus, according to the fourth moments, 
we only know that the J„{bys (0 < 6 < i) are different and the 7j^(6)'s 
(i < 6 < 1) are different. 

An added note. After the paper as submitted we learned from the As- 
sociate Editor and the referee about a related preprint "Limiting Spec- 
tral Distribution of Some Band Matrices" by Basak and Bose at the site 



http://www.isical.ac.in/~statmath/html/publication/techreport.html Al- 



though their paper and ours contain the same results for Hankel and Toeplitz 
band matrices, the former assumes less integr ability on the entries of a ma- 
trix, allows more general "rates" for the bandwidth, and also covers more 
ensembles of "structured matrices" related to Toeplitz matrices. On the 
other hand, our paper covers Hermitian Toeplitz matrices and gives low or- 
der moment calculations. Besides, our paper gives a different method for 
analyzing Toeplitz matrices by treating them as a linear combination of de- 
terministic matrices with independent coefficients. This method can be used 
to derive other results, including those that deal with semicircle law. 
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